CONSECUTIVE MINORS FOR DYSON'S BROWNIAN 

MOTIONS 

MARK ADLER, ERIC NORDENSTAM, AND PIERRE VAN MOERBEKE 

Abstract. In 1962, Dyson |Dys62| introduced dynamics in random matrix 
models, in particular into GUE (also for /? = 1 and 4), by letting the entries 

^^ ■ evolve according to independent Ornstein-Uhlenbeck processes. Dyson shows 

the spectral points of the matrix evolve according to non-intersecting Brownian 

p^ ■ motions. The present paper shows that the interlacing spectra of two consecu- 

,„.^. five principal minors form a Markov process (diffusion) as well. This diffusion 

^S^ consists of two sets of Dyson non-intersecting Brownian motions, with a specific 

interaction respecting the interlacing. This is revealed in the form of the gener- 
ator, the transition probability and the invariant measure, which are provided 

OO ■ here; this is done in all cases: /? = 1, 2, 4. It is also shown that the spectra of 

three consecutive minors ceases to be Markovian for /3 = 2, 4. 

^ ; 1. Introduction 

In 1962, Dyson |Dys62| introduced dynamics in random matrix models, in par- 
ticular into GUE, by letting the entries evolve according to independent Ornstein- 
Uhlenbeck processes. According to Dyson, the spectral points of the matrix evolve 
>- 1 according to non-intersecting Brownian motions. The present paper addresses the 

^ ■ question whether taking two consecutive principal minors leads to a diffusion on 

^ . the two interlacing spectra of the minors, taken together. This is so! The dif- 

O i fusion is given by the Dyson diffusion for each of the spectra, augmented with a 

C^ I strong coupling term, which is responsible for a very specific interaction between 

the two sets of spectral points, to be explained in this paper. However the motion 
induced on the spectra of three consecutive minors is non-Markovian, for generic 
initial conditions. A further question: is the motion of two interlacing spectra a 
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determinantal process? We believe this is not the case; but determinantal pro- 
cesses appear upon looking at a different space-time directions. These issues are 
addressed in another paper by the authors. 

During the last few years, the question of interlacing spectra for GUE-minors 
have come up in many different contexts. In a recent paper, Johansson and Nor- 
denstam |JN06] . based on domino tilings results of Johansson |Joh05aj . show that 
domino tilings of aztec diamonds provide a good discrete model for the consecutive 
eigenvalues of GUE-minors. In an effort to put some dynamics in the domino tiling 
model, Nordenstam |NorlOj then shows that the shuffling algorithm for domino 
tilings is a discrete version of an interlacing of two Dyson Brownian motions, in- 
troduced and investigated by Jon Warren [War 07] . One might have suspected 
that the Warren process would coincide with the diffusion on the spectra of two 
consecutive principal minors. They are different! 

Non-intersecting paths and interlaced processes (random walks and continuous 
processes) have been investigated by several authors in many different interesting 
directions; see e.g. JNF98], jJ^An2].|John5a], [JohOSbj . |TWn4] . |KTn4] . jCTGna] . 
|MUW09], pefOSaj, [KT04] . [Kg09] . |X^^M05] . just to name a few. In particu- 
lar, in |TW04[ lAvMOSj . partial differential equations were derived for the Dyson 
process and related processes. 

The plan of this paper is the following. We state precisely all the results in 
Section |2j Some usefull matrix equalities are derived in Section [3] which are used 
in Section H] to derive transition densities for the various processes considered. 
Stochastic differential equations are derived in Sections [5] and El The fact that 
the the spectra of three consecutive minors are not Markovian for generic initial 
conditions is demonstrated the last Section. 

There is a companion paper by the same authors aiming at determining the 
kernel for the point process related to the Dyson Brownian minor process along 
space-like paths |ANvM10] . 

Acknowledgement: PvM thanks S.R.S. Varadhan for several insightful conver- 
sations in the beginning of this project. 



2. The Ornstein-Uhlenbeck process and Dyson's Brownian motion 

Consider the space Tin oi n x n matrices B, with entries Bki G M, C, EI 
(/3 = 1, 2, 4) satisfying the symmetry conditions 

(2.1) Bm = i?4- 

Any element 2; G M, C, EI admits a decomposition z = z^^^ + X]r=i ^^^^^r, with e^'s 
satisfying 

el = el = el = -1, 6162 = -6261 = 63, 6163 = -6361 = -62, 6263 = -6263 = Ci. 
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The conjugate * of an element z G M, C, EI and its norm are given by 



/3-1 



z = z 



(0) 



/ ^ ^ ^Tl 



\z\ = zz 



E 



Ar)2 



r=l r=0 

with z admitting a polar decomposition z = \z\u, with |mP = X]r=o ^^^'^'^ ~ -'-• '^^^ 
matrices B G "Hn , as in fl2.ip . correspond to: 

{real symmetric n x n matrices, for /3 = 1 
complex Hermitian n x n matrices, for (3 = 2 
self-dual Hermitian n x n "quaternionic" matrices, for (3 = A 



with the compact groups of vol(Wn 



{/3)^ 



Unitar | (2.2) 



0(n), P = 1 

Ui^^ := lu{n), (5 = 2 

Symp(n), /3 = 4, 



acting on it by conjugation. 

For /3 = 4, it is well known that the quaternionic entries z can be represented 
as follows 



quat (2.3) 



z = z 



(0) 



/3-1 

+ E 

r=l 






2^°) + iz 



(1) 



z(2) + 2Z 



(3) 



.^(2) +,^(3) ^(0)_,^{1) 



So, the nxn quaternionic matrices B can be turned into 2nx2n self-dual Hermitian 
matrices 5, of which the real spectrum is doubly degenerate. Here, we shall define 
the n distinct eigenvalues as the spectrum of B. Unless stated otherwise we shall 
be working with the nxn quaternionic matrices, rather than the 2n x 2n Hermitian 
matrices. Also, when working with matrices having quaternionic entries, the trace 
will be defined in the usual way, that is as the sum of the diagonal entries of the 
n X ri-matrix. 

The determinant of a matrix B G "Hn is given in terms of the 2n x 2n matrix B 
(as defined in (12. 3p ). by the following procedure: first define the skew-symmetric 
2n X 2n matrix B by the following product: 



bold-b 



(2.4) 



B- 







-I 



®In 



2.24' 



and then "det 5" is defined as 

i 
(2.5) det5 := Pfaff(B) = {dei{M)Y'^ = J^i-'^T'^YlB^^Bp., . . . Bs 
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where p is any permutation of the indices (1,2, . . . ,n) consisting ot i exclusive 
cycles of the form (a — )■ /3 — )■ 7 — )■ • ■ ■ — )■ a); see Mehta (Meh88j . In particular, 
this means that 



quatdet (2.6) det(A/ - 5) = ]^(A - Aj), spec 5 = {Ai, . . . , Ap}, 



with the Aj being the double eigenvalues of B. 

The following normalization constant Z~\ will come back over and over again: 

(2.7) Z„J:=2-t(^^j , withiV:=iV„,^:=- + ^n(n-l). 

Dyson's idea was to let the free parameters of the matrix evolve according to 
the SDE (Dyson process) 

[2 
dBii = -Biidt + J -dbii, z = 1, . . . , n 

V r^ 



1.3| (2.8) 



dB^f = -B^fdt + -^dblf, i, J = 1, . . . , n and £ = 0, . . . , /3 - 1, 



where dba and db\j , for I < i, j < n and / = 0, . . . , /3— 1, are independent, standard 
Brownian motions. Since the Ornstein-Uhlenbeck diffusions are independent, the 
Dyson process on the matrix B has a generator, which is just the sum of the 
OU-processes above: 



/I 9^ d \ / 1 (9^ 

m^ (2^9) .Ao,, := 5: (^-^ - B,— j + Y. 5: a 






OU-tr 



with transition probability, setting c := e * and using the constant (12. 7p . 

P[Ei G ci5 I Bo = 5] =: pit, B, B) dB 
(2-10) Z-l, 



«./3 .-^7T^Tr(B-cB)2 



'1 _ c^)N„,i3 



e-^(T^^ ^'^"-'-"^ t/S, 



where rfi? is the product measure over all the independent parameters Ba, B- . 
The transition probability (I2.10p satisfies the Fokker-Planck equation 

(2.11) ^=^DysP, 



with 



FPl 



2/iA<9,«5i ix-^^a,«a 1 



pi2) Ai,.^U\Y.4rk'~,A E E^"^ 



\ 4 = 1 1<«<J<" t=0 IJ 13 
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with a delta-function initial condition, p(t, B, B)\ = 6{B, 5), and with invariant 
measure (density) 



Invl 



(2.13) lim p{t, B, B) = Z-\{h{B)Y , with h := h{B) := e 



-iTr B2 



i— >oo 



Dyson discovered in |Dys62| the surprising fact that the process restricted to 
spec(-B) := {Ai, A2, • • • , A„} is Markovian as well. This is the content of Dyson's 
celebrated Theorem (Theorem 12.11) . 

Before stating the main Theorem, we define diagonal matrices X = diag(xi, . . . ,Xr, 
and Y = diag(|/i, . . . , ?/„), vectors w,v E M", C", H" and the inner-product 
{w , v) = Xli ^j'^i*- Then consider the integral 



1.21 



(2.14) G^^{X,Y-w,v) := / ^^ rff/e(Tv^t/yc/-i+2ReKc/.))^ 

and its integrand 
gI^ (2.15) e^f)(f/;X,r;^,t;):=e(T^^^^^''+2^^<"''^^>). 

For w = V = 0, this is the more familiar integral 

(2.16) Fi^\X,Y) := Cf (X,F;0,0) = / df/e^--^^^^"' , 

which for /3 = 2 gives the Harris-Chandra-Itzykson-Zuber formula: 



HCIZ 



^9 17^ rC2)(x Y\-^^^^^—\l^hi^'n 



Tl . 



Does the integral (I2.14p admit such a representation ? This is an open problem. 

In the following Theorem, formulae (I2.18p . (I2.20p and (I2.22p are due to Dyson 
|Dyi62| . 



ThrDyson Theorem 2.1. The Dyson process restricted to its spectrum spec{B) = A := {Ai, 
. . . , A„} is Markovian with SDE given by: 



1.4 



h^s^) 



(2.18) d\i= {-\i + }_^^-^\dt+J^dbu, i = l,...,n. 



1.20 



1.5' 



3.13" 



footC 
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Its transition probabilityj, with c := e~*, 



(2.19) P[Ai G d\\\o = A] =p^{t,\,\)d\i---d\n 



^-1 



(1 _ c2)^n.^ 



5^E?(A?+c2A?)^(^) 



/3c 



A,A)|A„(A)|'^n^^i 



satisfies the Dyson diffusion equation, with delta-function initial condition (px\ „ 
6{\, X)) (forward equation) 



(2.20) 



dp^ 



dt 
The generator is 



-4>.. ^^hA::=l±^m^)f^^j^,- 



i=l 



(2.21) 



A, 



El-— + 



/3 3Af 



-^'+Ev^ 



Aj — Aj 



9A" 



and the invariant measure of the Dyson process on B, projected onto spec (5), is 
given by the GOE(n), GUE(n), GSE(n) measure for P = 1, 2, 4 respectively: 



lT28] (2.22) C-U^niX)fdX^---dX^, «;zi/i$„(A)=e-5Er^,|A„(A)|. 



For completeness we shall prove (^A^ . fl2:22|) in Section land fl238D . fl2:20|l in 
Section [51 

It is remarkable that the Dyson process is not only Markovian upon restriction 
to the spectrum of any single principal minor B, B^'^~^\ i?^""^), . . . , of sizes n, 
n — 1, n — 2, . . . , but also upon restriction to any two consecutive principal minors, 
in particular, 

(spec B, spec 5*^"^^^) := (A, /i) := ((Ai, . . . , A^), (/ii, . . . , /in-i)), 

with intertwining property 

(2.23) Ai < /ii < A2 < /i2 < ■ ■ ■ < A„_i < iin-i < Xn- 

We denote by Ax and Afj, the generators of the consecutive spectra spec B and 
spec B^"'~^\ as defined in O2.20p . Define the characteristic polynomials of the two 



C^ n is the norming constant for the Gaussian ensemble for general /3, as obtained from the 
Selberg formula (see Mehta |Meh04] . formula (3.3.10)), (see (E?]) for N) 

rfi + 



C-^ = (2vr)-t/?^n 



j=i 



r 1 
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consecutive minors B and B^'"'~^\ 



TIO] (2.24) Pn{x) = l[{x-Xa), 



a=l 



n-l 
Pn-lix) = Y[ix- fip), 
/3=1 



vdm 



and the Vandermonde determinants 

An(A) := n(^.- - ^*) > 0, 

j>i 

yZ,.Z,0) n n—1 n 

A„(A,/i) := nn(^^ --"i) = IlPn-l{\) = l[Pn{^i), 



n-l 



i=l j=l 



n{n—l) 



conjl 



with A„(A,/i)(— 1) 2 > because of the intertwining. 

In order to state Theorem 12. 2[ we need the following property of any matrix 
B e T-Ln , not only can B be conjugated by a matrix f/*^") G Un (see (12. 2p ). in 
the standard way, such that 

(2.26) ([/("))-^5f/(") = diag(Ai, . . . , A„), 



but also by a matrix of the form 







^J, with t/("-i) G Wj^^, to yield 



2.22 



bordered matrix -Bbord '■ 










(2.27) 












/ /"i 


•• 





nui \ 


V V V V ~ 


u 


/i2 •• 





r2U2 




u 


■■ 


/^n-l 


Tn-lUn-l 




ViM* 


r2U*2 ■ ■ 


■ r^-i<~i 


rn J 



--■.B 



bord 



with \ui\ = 1 (angular variables) and with r^ > for 1 < i < n — 1 and r„, given 
by 



S (2.28) 



PnjfJ'k) 
PLlM 



n-l 



> 0, 1 < A; <n-l, 



'•n :=^Ai-^/ii. 



The conjugation in (I2.27P transforms the last column f of i? in the last column of 
the bordered matrix -Bbord (except for the last entry); i.e., 
E (2.29) 

Jj{n 1)^ ^ {nUi, ..., Vn-lUn-lV, and Bnn = rn, with V := {Bi^n, ■■■, Bn-l,nV. 

These facts, (I2.27p . (I2.28P and (I2.29p . will be discussed and shown in Section [31 

The next statement is the analogue of Theorem 12.11 for the case of the spectra 
of two consecutive minors. 
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Th:Dyson2 Theorem 2.2. The Dyson process on B restricted to 

(spec 5, spec 5^""^)) = (A,/i) := ((Ai, . . . , A„), (/ii, . . . ,/i„_i)) 
is a diffusion (A(t),yu(t)) as well, with the following SDE: 



SDE-lambdaO] (2.30) 
d\ 



2 Pn_l(Ac 



-(-^"+Ei:^)*+yi^ 



{\ 



nV^ci 



E a/2 rjrj dbjj ^ rfdbu ^ a/2 rjdbjn ,, 



( 1^/2" 

\ ^^^ A^T ~ /^e / V P 

m terms of independent standard Brownian motions {dbii,dbij}i<i^j<n- Its tran- 
sition probability^ is given by: 

PXf, (t, (A,yu),(A,/i)) t/Ad/i 

=¥[{\t,^lt) e {d\,dfi) I (Ao,/io) = (A,/i)] 



1-i 



":/3 „-^7T^Ei(A?+c2A?) 



1.24 



(2.31) 



1-c 



,2W 



e 2(i-c^) ■ 



n-l 



(5/3-l)2(n-l) 



l[dfl^^-'\u,)dQ^^~'^ 



.Ui 






I3c 



rm. 



\n-l 



— C 



2 V «"-«;! ; V I'^iJl 



rm, 



\n-l 



l^crnTn 



X 



1.12' 



n n—l 

e^^T^ I A„(A) A„_i(/i) I I A„(A, /x) |(^i) J] rfA, J] d^i,, 

1 1 

where the ri 's are given by 112. 28\} . It is also a solution of the following forward 
diffusion equation, with delta-function initial condition 

dpx^ 



1.12" 



(2.32) 
where 

(2.33) 



'-Am' 



^I,: 



^ A PA/., with A' ■.= Al+A'^+Al 

1 Pn-l{\) Pnifij) 



2 A^ 9 d 



/StrS^^^^-^AlA.-/",)^ ^^(A.) P^i(/i,) 



f ootZhat 



The constant reads 



^-1 _ 



/3^-''(r(i + §))"-! 



(27r) 2 ^l("-i) n;ji' r(l + f )(vol(5^-i))"-i 
and vol(5'=) = STrHtTi^ (2 /J'''^ (cos 6l)M6l) for k > 2, vo\{S") = 1 and vo\{S^) ^ 27r, which is 



proved by induction on fc; so vol(S' ) = 47r and vol(5 ) = 27r 
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1.29 



restr 



general 



Cor :unis 



and where A^ and ^^ are defined by f l2.20p . The Dyson process restricted to (A, /i) 

has invariant measure, (see /i2.25]) ). 

(2.34) 



n-l 



1 1 

The SDE fl2.30p and generator (12.321) are computed in Section [6] while the espres- 
sions for transition density (I2.3ip and invariant measure (I2.34p are proved in Sec- 
tion H 

Note that it is an immediate consequence of Theorem 12.11 that the generator 
^Dys; defined in (12.90 . acting on the Aj and /ij, has the form 

(2.35) ^Dys(Ai) = A(Ai) and ABysifJ'i) = A^{fii), 

where ^a and A^ are defined by (12.211) . 

Whereas all statements in this paper hold for P = 1,2,4:, a, part of it can be 
extended to general /3 > 0, as will be shown in section |6l after the proof of Theorem 
E21 



Corollary 2.3. For general /3 > 0, the SDE 112. 30\) . in terms of the independent 
standard Brownian motions {dbu, dbij}i<i<:j<n, defines a diffusion, whose generator 
is given by the same equations l[2. 3^) . and whose invariant measure is given by 
l[2.34\) - Moreover, this diffusion restricted to the Aj 's (or to the /ij 's) is the standard 
Dyson Brownian motion i\2.18\} . 

The following corollary shows that the /ij's in Aj < /ij < Aj+i are repelled by the 
boundary and fluctuate in unison with the boundary points, when they get close. 



Corollary 2.4. The nonnegative gaps /ij 
satisfy, in the notation of Ii2.30\) . 



Aj and Aj+i — ^i for 1 < i < n 



unis 



(2.36) 

with 
(2.37) 



d{fii - Xi) = Fi{\, n)dt + ^J y,i - Aj ^ akidbki 

l<k<e<n 

d{Xi+i - fii) = Fi{X, ^)dt + v^Ai+i - /ii ^ akidbke 

l<k<e<n 

some aki = 0{1) for /ij ~ Aj and some dij = 0{1) for jii ~ Aj+i. 



F^{X,^i)\ 



>0, F,(A,/i)| 



>0 



This is to be compared with the Warren process [War 07] . which also describes 
two intertwined Dyson processes A and /i, but with an entirely different interac- 
tion: namely the /Xj's near the boundaries of the intervals [Aj, Aj+i] behave like the 
absolute value of one-dimensional Brownian motion near the origin. 
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As we saw, the Dyson process on B, restricted to the spectrum of one principal 
minor or the spectra of two consecutive minors leads to two Markov processes; 
opposed to that, we have the following statement, which will be proved in Section [71 



Theo:3minors| Theorem 2.5. The restriction of the Dyson process restricted to the following 
data 

(spec E, spec 5("-^\ spec fi("-2)) := (A,/i,z/) 

is not Markovian for generic initial conditions on B, i.e., the joint spectra of any 
three neighbouring set of minors of B are not Markovian, for (3 = 2 and 4. 

3. Some Matrix Identities 



sec: matrix 



2.16 



2.28 



2.24 



2.25 



In this section, we prove formula (12.271) for r^ as in (12. 28 p . In the course of doing 
that we will also prove the following formulas: 

71—1 2 1 / " n— 1 \ n— 1 

(3-1) E''? + t = ME^?-E''? 

1 \ 1 1 

One also has the (often used) identities 



and TT 



|An(A,/i)| 



n-l 



(3.2) Y^ —^ — + r„ - A£ = and 



Ki^t 



/i. 



n-l 



(A, 



n-l 

E 



2. 15 I (3.3) 



j=l i=l 

Finally, one has, for fixed (/ii, . . . , /i„_i) and fixed (mi, . . . , m„_i 

„-i A„(A) " 



A^ — /ij 



+ 1. 



\{dr]drn = {-iy 



\{d\i 



1 A„_i(/i) ^ 

Proof. From the form of the matrix -Bbord as in (12.271) . one checks (see (I2.25P and 
also the formula (12 .Sp for the determinant in the quaternionic case) 

n 71—1 /n— 1 



(3.4) J](A, - A) = det(5bord - A/) = \{{^i^ - \) [Y, y^ + r 



.i=l 



A 



/Ltj 



from which it follows that □ 



(3.5) 



n-l 

E 

j=i 



rp^ 



A -/ij 



r-n + A 



Pn{\) 
Pn-liX) 



= X - (ai(A) - a,{f,)) - J (ai(A)o-i(/.) + aM - a,{X) - aUf,)) + O(^). 

Then taking residues in formula (13. 5p yields the first formulae (12.281) and thus the 
formula for Yli "^f i^ (13. ip . Comparing the coefficients of A° and the A~^ on both 
sides of (13 -Sp yields the first formula of (13. ip . Setting A = A^ in the expression 



the o'fc(A) are symmetric polynomials: cri(A) — ^^ A^, (T2(A) = X]i<i -^i-^ji 6tc. . 
for (7k{p) 



The same 
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2.27 



2.29 



(I3.5p and its derivative with regard to A implies the two sets of n identities f l3.2p . 
in view of the definition of P„. Formula fl3.3p amounts to computing the Jacobian 
determinant of the transformation from Ai, . . . , A„ to ri, . . . , r„; to do this, take 
the differential of the first of the n expressions appearing in (13 ■2p (as functions 
of Ai, . . . , A„ and ri, . . . , r^), keeping the /Xj's fixed and use the second of the 
expressions (13. 2p : 







(3.6) 



n-l 

E 

1=1 

n-l 

E 

i=l 



dr'j 



n-l 



Af — /ij 
dr'j 



Xf 



^J'i 



dr„ 



+ drn 



i + E 



^ {Xe - fiiY 



dXf 



Pn-liX) 



dXf 



which in matrix form reads 
/ drf \ 

r 



\ dVn J 



diag 



PLi^i 



PLi^ 



nV^n) 



n-ll'^l 



(Ai 



-Pn-l(A^ 



/ rfAi \ 
dX2 



\ dXn I 



where (by Cauchy's determinantal formula) 
(3.7) 

/ 1 1 ... 1 l\ 

Ai— /ii Ai — U9 Ai — /i^ — 1 



1 



r := 



Al-/t2 
A2— Ml A2— /12 



\A„— /^i A„— /i2 



Al— /in-1 
A2— Mn-1 



with det r 



n fJ'Ti — l / 



^Yn-l)(^ + l) An(A)An_i(/x) 
A„(A,yU) 



An~M' 



The formula (13. 7p for the determinant follows from the observation that det F has 
homogeneous degree 1 — n and vanishes when A„(A)A„_i(/i) does and blows up 
(simply) when and only when A„(A,/i) vanishes. Thus we have 



d{rl 






n#^(detF)-i = (-l)"-i ^"^^^ 



a(Ai,...,A„) -'-ip^_^(_x.)^^""^^ ^ ^^ A„_i(;u)' 

This concludes the proof of formulas stated in the beginning of this section. 



D 



transition 



4. Transition Probabilities 



A quick review of the Ornstein-Uhlenbeck process (see Feller |Fel71] ) : it is a 
diffusion on M, given by the one-dimensional SDK, 



pnl (4.1) 



dx = —px dt -\ — -= db, 

V/3 



12 MARK ABLER, ERIC NORDENSTAM, AND PIERRE VAN MOERBEKE 

and it has transition probability (c := e"''*) 

P[xf & dx \ Xq = x] =: povit] X, x) dx 



p/3 



7l(l — C 



1/2 
,2N / e^P 



p/3{x — ex 

1-C2 



;^^2 



(ix. 



The transition probabihty is a solution of the forward (diffusion) equation, with 
5-function initial condition u 



T3] (4.2) 



<9pou 



1 d^ 



d 



dt V 2/3 (9x2 Q^ ^ 
and invariant measure (density) 



px) Pou 



a 



a 



2/3 V5x'^'"'^''^ax(/)^^(x) 



Pou, 



T4] (4.3) 



Via^j 




-p/3a;2 



lim j9ou(i;a;,x). 



T5 



Proof of transition probabilities l[2.1(J\) . \2.19^) and \2.31]) . 

(i) The Fokker-Planck equation for the transition probability of the Dyson process. 
The Dyson process consists of running the free parameters of the matrix B e 
"Hn , as in (12. ip . according to independent Ornstein-Uhlenbeck processes, with 
p = 1, the diagonal with (3 — )■ (3/2 and the off-diagonal parameters with f3 ^ (3. 
Remembering the definition (12.71) of A^ = Nn^i3 and the definition of the trace (after 
(12.31) ). one has, setting c = e~*, and using f|2.8p . (14. ip and (14. 2p . the transition 
probability for the Dyson process is given bjo 



(4.4) 



/3-1 



j=l 



p{t,B,B)=l[pov{t;B,„B,,) J] l[povit] b';^ ,Blf) 

l<i<j<n e=0 

n((2vr(i-c2)//3)i),<n^„n 



I a 'J U- 



V 



>(l-c2)//3)5 



\ 



1 



Af„, 



2"/2(|(l-c2))''"' 



-e 2(1-='*) 



Tr(B-cB)2 






'1 _ c2)^",/3 



H^Tr(B-cB)2 

-g 2(l-c2) '^ ' 



The backward equation becomes the heat equation with {x,t) H> (xe''*, ^-^ — ) 



with constant Z. 



n,/3 



2-«/2(|)-A^„,^. 
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yielding fl2.10p . while \iiat^ooP{t, B, B) = Z^^(/i(i?))^ is immediate, showing 
fl2.13p . Moreover, from (14. 2p . one computes for p{t] B, B), 

(4,5) 



|p«;5.B)=| 



. n 



i=l 



d 



2\Bii^ 



/3-1 



d 



1 



OB, 



a VB 2\J->ii 



Bi, 



+ 



z E E 



d 



d 



4 ^ ^ ^n^^i'^^^^'i'' ^ni^) ^jni^h 






dB-:> MBl^' 



p{t;B,B) 



(4.6) 



''lj:^m-' ' 



(3 



2 ^ 95, 



dBu KB) 



/3-1 



+ 4 2^ Z^^rW 



« h(B) " 






p{t;B,B) 



with 
(4.7) 



h{B) = constant x ]^0_g/2(-Bn) Y\_ 



ij ) 



i=l 



0<f.<P-l 



proving (I2.12p . 

(ii) The transition probability Ii2.19\) and the invariant measure for the Xt process. 
Set 

A = (Ai,...,A„) and A = (Ai, . . . , A„). 

By the Weyl integration formula, given B = U\U~^ and initial condition B = 
UXU^^, express dB = d{UXU^^) in formula (I2.10p in terms of spectral and an- 
gular variables dB = Zn,i3C~^p\An{X)\^dU YYi dXi, with Haar measure dU on Un 

normalized such that Vol(Wn ) = 1, with C~^ defined in footnote [H This yields, 
using the transition probability (14. 4p . 



2T8] (4.8) F[BtedB\Bo = B]=p{t;B,B)dB 



^-1 



[1 - C2) 



,2W 



g-,(T^Er(A?+c^A?) 



X ei- 



■TrU\U-^U\U-^ 



AMfdUHdX,. 



Note that the constant C~\ is compatible with the fact that for t — )■ 00 this 
transition probability tends to the GUE-probability; see below. 

We now compute the transition probability px(t; A, X)dX for the spectrum of the 
Dyson process; this will be a model to compute the transition probability for the 
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(A(,/ii)-process. So, consider 
(4.9) 

Px{t; A, X)dX = P(Ai G dX\Xo = A) 



/ ,,, /_ ,,, F[B,EdB\B, = B]F[B,EdB\speciB,) = X] 



using fl4.8p above, using the following conditional probability formula: 

r - 1 . s -1 FlB.edB, spec(fin) e dX] 

(4.10) P i?o G di? spec(i?o) = A = ^ ° . ,^ M ,°-,i ' = dU 

P [spec(iJo) G dX\ 

and finally using the integration (I2.14p . 



thus yielding fl2A9|) . 

Letting t — )■ cxd, (equivalently c — > 0) in (14. 8 p proves formula (I2.22p for the 
invariant measure, taking into account that Fn (0, F) = vol(W^) = 1. 

(ii) Proof of the transition probability 1(2. 3 1\) and the invariant measure for the 
{Xt, fit) process. The proof of (I2.3ip in Theorem 12.21 proceed along similar lines. 
First observe the identity 
condP I (4.11) 

F[Bt G dB I (spec(5S")),spec(i?i"-'))) = (A,/i)] 
= /_ P [Bt EdB \ Bo = B] P [^0 = d5|(spec(5j"^),spec(5j""^^) = (A,/i) 

with 

nl^\X,fc) = 7/(f) n {(spec(i?W),spec(i?("-i)) = (A,/i)} 

Next we compute the two probabilities in the integrand of the integral (14. lip : 
(i) The first integral equals p{t; B, B)dB, as in (14. 4p . Since the Haar measure dB 
is the product measure over all the free parameters, one will express dB as the 
product of Haar measure dB^"'~^^ on the (tt, — 1) x (n — 1) minor and the measure 
n i<i<"-i dBlJdBnn on the last row and column, remembering the expression fl2.7p 

o<e</3-i 



condl 
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for Nn^p, thus giving: 

(4.12) 

nSt e dB\B, = B]= . ^"\^.^ y ^i^^'^^'-'^'^^'dB 

i — c ) "•f' 






n-1/3-1 
i=l £=0 

As mentioned prior to (14. 8p . one can set in (I4.12p . 

71—1 



0<t<l3-l 



Bo] (4.13) dS("-i) =Z„_i,^C-ii,^|A„_i(^)|'^df/("-i) J]rf/i,. 



i=l 



Bl 



In (I2.27p . it was shown that upon conjugation by an appropriate matrix f/^" ^^ G 
W^_\, the matrix B could be transformed in the bordered matrix -Bbord, as in (I2.27P 
and (I2.29p . with (riUi, . . . ,rn-iUn-i)~^ = U^"'~^'^v and \ui\ = 1. Using the same 
inner-product as in the formula just preceding (12.140 . but for n — 1- vectors, and 
using the associated norm || ||, one finds, using the above, 
(4,14) 



l<i<n-l 
0<f</3-l 



If — cvW^ 



||;7("-i)^||2 + cl?7("-i)f II' - 2cRe([/("-^)t;, U^^'-^^v) 

||;7(""i)^;||2 + c^\\u^--^^vf - 2cRe(t/("-^)f , ([/("-^)([/(""i))-^)(f/("-^))f)) 

n— 1 n— 1 

JZ^' + c'J^r- 
1 1 

- 2cRe((riMi, . . . , r^^^Un-iV , {U^''-'\u'~''-^^)-^){nu^, ..., f„_iM„_i)^) 
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Given that f/("~^) is fixed and that det{U'^^^^^) = 1, and since the expressions Ui in 
fl2.27p have \ui\ = 1, the differential below can be written in terms of a product of 
differentials d{riUi) along S^^^ C M'^, expressed in polar coordinates, thus yielding 
differentials involving the r^'s and volume elements on the unit sphere 5*^"^: 



n-l/3-l 



n— 1 n— 1 



B2 



n n ^^in = n ^^^ = n diiu'^^'-'^r^nu,, ..., rn-mn-iV), 

i=l i=0 1=1 i=l 

n-1 

(4.15) =|det(f/("-^))-^in^M.) 

1 

n-l 

= TT rf ~ drid^l (ui) and dB^n = dr^. 



coiid2 



Thus all together, setting 005]) . JKT^ and IKT^ in ( CT^ . we have shown that 
(4.16) 

^-1 



e .(T^EE^(A«?+-^M?) 



(1 _ c2)A'n-i,/3 



n-l 



X 



X 



^ni^n-1,/3 _l^l(j2r'rHhi)+-HU-'rj^fl)) 



(1 _ (.2-)(Af„,^-Af„_i,^) 



e (1-=^) 



X e 



^Re((r,«,)?^-i,C/("-i)(C7("-l))-i(f,S,)5'-l> -^r„f, 



n-l 



gTZ^'-n^-n ^^^ -Q ^/3-l ^^. ^^ 



(/3-1), 



M,: 



^""^^"-^'^ e-^(T^^"^^?-^^^\T^-1A„(A)A„_,(;.)| |A„_,(A,/.)|f- 



(l-c2)^-^a-l,/3 

X ^f_\ff/("-^)([/("-i)) 



(/3) f TT(n-l)^rT(n-l)\-l. P^ ,. -.. P'^ ^„ „. \n-l ^^ „- Nn-1 



r/^>/^; 



X 



rf^(n-i) Y[da 



n— 1 n— 1 n 

(/3-1) 



-(riMiji ,l^i^**;i 



Uii 



Y\ dfii Y\ d^i 



1 



1 1 



In the last equality we have used identities (13. ip and (13. 3p and the definition f l2.15p 
(ii) The second probability in (14. lip takes on the following value: 



cond2' 



(4.17) 



P 



(n) 



An) 



B^""' e ciE (spec(5^"0,spec(5^""'0 = (Kf^) 



(n-l)^ 



dU^""-^^ 



(vol(5/5-i))«-i 
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Indeed, the probability fl4.16p . when t — )■ oo ( which amounts to letting c — )■ 0), 
tends to the invariant measure for the Dyson Brownian motion; instead of the 
usual representation in the variables {\i,U^"^), this gives the expression of the 
GUE-probability in the variables {Xi, fj.j,Uk,U^^~^^) : 



P[5(") e dB] 



Z-\e~>^' 



ri,;S 



dB 



?(") 



(n) 



S (4.18) 



lim ¥[Br e dB\B^''' = B 



t— ^oo 



n-1 



n-1 



fnA>^, /i)rff/("-i) n "^^^''^ ("') n ^/^^ n ^^^ 



'(/3) 



^ „-i \n-l\ 



with (using ^;,1\(f/; 0,/i; 0, (f.^O? 



(4.19) /n,/3(A,/i) := %^^^^e-f ^"^?|A„(A)A„_i(/.)| |A„_i(A,/.)|^^ 

C-n-1,/3 



, n-1 



This also shows that 'Hi^\\, fx) ~ U^^^ x (S^'^-^))" '. Using ( HTTSj) . one checks 

that the conditional probability equals 

(4.20) 



P 



5i") G rf5|(spec(5^"0,spec(5^""'0 = (Kfi) 



P 



?(") 



?(")^ spec(5i"-')) 



5^"^ G dfi, (spec(5^"0,spec(5^"~'0 e (^A, rf/i) 



P 



(spec(5j"^),spec(5i""^^) G (rfA, d/x) 

/„,^(A, /i)rff/("-i) nr' dn\^''^ (w.) nr' ^/^^ n" ^a. 



/n,/3(A,/i)ni rf/^iHlC^A, 



wf_\x(s(/'-i)) 



_^d^(«-i) nr'^^f"'^' 



u,: 



(vol(5(/3-i)) 



n-1 
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3.10' 



3.12 



confirming expression fl4.17p . Setting fl4.16p and f l4.17p in f l4.1ip and using the 
integral f l2.14p and the identification just after (I4.19p . one computes: 
(4.21) 
Px,f,{t; {X, fi), {X, n))dXdij 



P 



Btenif\x,fi) JB&U^^Wfi) 



An) 



>H 



B',> e dB B\;'' = B 



P 



^S") e d5|(spec(5S")),spec(5f"~'^ 



;-(^./^) 



1^ V"~^l'\2 I r'2\2\ ^c 






1 — C^ 1 — C^ 

n-1 



n-1 



n-1 



iA„(A)A„_i(/i)i iA„_i(A,/i)i^i n^^*^''''^^"*) n^^s''~'^(^') n^/^^n^^^ 



i'-i 



(T^E..i (K+c ^JeTr^^^^"|A„(A)A„_i(/i)| |A„_i(A,/i)|i-^rf/irfA 






(3c 



n-1 



-(r^Mi)^ \(riMi^ 



)rM l[dn[^-'\u,)da 



^'-'\u. 



were we used the translation invariance of dU'^'^ ^^ and vol(W^_\) = 1; also one 
checks the value of the constant Z~\ to be the one given in footnote [2j This 
establishes formula fl2.3ip for the transition probability of the (At,/if)-process. 

The statements concerning the invariant measures, f l2.22p and f l2.34p follow im- 
mediately from f l2.19p . fl2.14p . fl2.3ip . by letting t — )■ cxd in the transition probabil- 
ity. This concludes the proof of the formulae for the transition probabilities fl2.10p . 
fl2.19p and invariant measure fl2.34p , appearing in Theorems 12.11 and 12.21 D 

Remark 4.1. The diffusion equation f l2.32p . which will be established in section [5l 
can also be used to confirm the form of the invariant measure, at least for /3 = 2. 
On general grounds, the density of the invariant measure, namely 

(4.22) h^.{X,ii) := Ce~^^'^^"\AniX)An-i{fi)\ |A„(A,/i)|i-\ 
is a null vector of the forward equation, i.e. 

(4.23) ^^/A,(A,/i) = {Al + Al + Al^)h,{X,fi) = 0, 
with A defined in (12.321) . For /3 = 2, more is true; namely 



"3T13] (4.24) 



a!{X)i, 



nin — 11 



A/i 



'A/i5 



^ZW. 



nin — V. 



Xfj, 



'■Xfi- 
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Once this is shown, it follows that 

(4.25) Aj^h,{X,fi) = -n{n - l)h,{X,fi). 



So it suffices to prove (14.241) . First observe that A„(A) and A„_i(/i) are harmonic 
functions, i.e. 



n / a \ 2 n—l / a \ 2 



and also homogeneous functions so that acted upon by the Euler operators. 



i:A.^A„(A) = ^^^A„(A), 



^ dX, 



n-l 



E9 ^ / X (^ - 1)(^ - 2) ^ , , 
/i,— A„_i(/i) = ^ 'j ^A„_i(/i). 



Now compute from (l2.2Up and (14.221) that (remember ^n{X) := e 2^1 -^i 

|A„(A)|) 



ATf.^r _ 1 V- -9 2 d A„(A)A„_i(^)e 

2 4-aA/ aA, "^^^ 



-E^A^^ 



ie-^^^=A„_,(,)|:(iJ-2A.|-)A„(A) 



n{n — 1) ^ 
n -'A/i 
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and also that 

d A„(A)A„_i(/i)e-^"^' 



1 """^ d 



dfii ($„_i(;U))2 



n— 1 rj 






y-n-l 2 d 

gZ^l A*! 



A„_i(/i)e" 



'Erv? 



|A„(A)e-Sl^^? (^g( 



92 



djjf d^. 



d 

2/ii-— )A„_i(/i) - 2(n-l)A„_i(/i) 



-A„(A)e-^"^?((n - l)(n - 2) + 2(n - l))A„_i(/i) 
n(n — 1) ^ 



This ends the proof of identities fl4.24p . 



secSEEd) 



3.1 



5. Ito's Lemma and Dyson's Theorem 

To fix notation we repeat some well known facts from stochastic calculus in a 
way that will be useful later. Given a diffusion Xt G M", given by the SDEu 

(5.1) dXt = a{Xt)dbt + a{Xt)dt, 

where dbt is a vector of independent standard Brownian motions, where x, a{x) G 
M" and (t{x) an n x n matrix. Then the generator of this diffusion is given by 

(5.2) -4=^E(^^'")^^(^)7)5:r + E«^(^) ^ 



«j 



dxidxj 



3T2] (5.3) 



and, by straight forward verification. 



dxi 



dX^dJ\j 



dt 



X 



The transition density p(t, x, x) is a solution of the forward equation (in x) 



3:3] (5.4) 



dp 
'dt 



A'p 



Moreover for a function ^f : M" i-)- W with gf G C^, the SDE for Yt = g{Xt) has the 
form 

3.51 (5.5) dn = J2^^dX, + lj2^dXdX,=J2iJ2^/^^)db, + h,dt, 

i i,j -' j i 



The subscript t in Xf and St will often be omitted. 
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for k = 1, . . . , p and for some function /i^; i.e., the local martingale part only 
depends on first derivatives of g. This follows from the standard multiplication 
rules of stochastic calculus [dtdt = 0, dtdb = and dbidbj = Sijdt): 



(5.6) 



dXi dXj = {tti dt + y^ an dhi) {aj dt + N^ ajk dhk) 

i=l k=l 

n n 

= (^criedbi)(^ajkdbk) = {(ya^)ijdt. 



(5.7) _ 

1=1 k=l 

More details can be found in any book on stochastic calculus, for example |McK05j 
or |0ksO3 . As a warm-up exercise, we first prove Dyson's original result, namely 
the formulae for the SDE and for the generator of Theorem 12.11 including some 
consequences. 

Proof of (I2.18P and (I2.20p in Theorem \2.1\ The Dyson process is invariant under 
conjugation by t/ G Uii ■, to be precise from (12. 9p . 

pit; UBU-\ UBU-^) = pit; B, B). 

Therefore, we are free, at any fixed choice of t, to reset 

Bit) ^ UB{t)U-\ for any U G U^^\ 

At any given time t, diagonalize the matrix B to yield diag(Ai, 
the perturbation 

diag(Ai,...,A„) + [dBy], 
where one defines the n x n matrix, for 1 < i < j < n, 





A„) and consider 



perturb (5.8) [dBi 



\ 



dBf-Y!t'dB 



V 



dBf+Y^UBf^e, 




and, for i 



with, by dlSD, 



[dBu] := diag(0. 



, CtJDii, 



w^T^in 



378] (5.9) dBi;> dBl'l = 6u'6j,>6u' 



dt 



'IJ — I'j' -Vf^Jf-Lt. n 



dBii dBjj 



,0), 



9A ^^ 



dBlfdBkk 



0. 



Remember by Ito's formula (15. 5p . one only needs to keep track of at most second 
order changes of the arguments. 
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Thus for non-diagonal perturbations (i 7^ j), one checkq^ 
= det fdiag(Ai, . . . , A„) + [dBij] - XI^ 



AM-Aa+C^A, 



aT'-i./Kcx 



"3T10] (5.10) 



JJ(A,-A) (A-A.)(A-A,)-5^(rfS 



W^2 






e=o 



Xh^Xa+dXa 



(non-zero function) x dXa, 

(non-zero function) x ((Aj — Xj)dXi — X]^=o^('^-^i 



/5-1mdWa2> 



for a^i, j, 
for a = i, 



(non-zero function) x ((A^ — Xi)dXj — X]£=o i^^ijY) ^^^ ^ ~ ii 

showing that an off-diagonal perturbation of the diagonal matrix B{t) = diag(Ai(t), . . . , A„(t)) 
yields 



dXi 



Etim. 



W\2 
ij ) 



Xi — A, 



-, dXj 



Xj — Xi 



and dXa = for a 7^ i,j. 



For diagonal perturbations {i = j), one finds 



3.13' 



detfdiag(Ai,. . . , Xn) + [dBii]~XI 



\t~^\a+dXa 



lliXe-X) (Xi + dBu - A) 



X^-^Xcx-^-d\, 



a-r^^a 



t^i 



(non-zero function) x rfA^, for a ^ i 

(non-zero function) x {dB^a — dXa), for a = i. 

and thus 

dXa = for a 7^ « and dXa = dB^a for a = i 
Then summing up all the perturbations, one finds the SDE (12.181) announced in 



Theorem 12. H using dB, 
dXi = dBj 



-B.:dt 



%dh^. 



—Xidt + \ r^dhii and formula (15. 9p . 






jV« 



Xi — Xj 



h^S^) 



dt+ \l-dbii, 



for i = 1 n. 



Then translating the SDE into the generator of the diffusion on (Ai, . . . , A„), one 
finds, by (15. 3p . that 



(5.11) 



A-i:(^J|+(-^.+i:^ 



)- 



Remember, for /3 — 4, the determinant is defined in 



SDE2 
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and thus 



Ai 



iV — ($ iX)f— ^ 




with $„(A) as in fl2.22p . confirming formula fl2.20p in Theorem 12.11 Finally ^Dys-^i = 
A^\i, mentioned in (12. 35 p . follows from the fact that the generator ^oys restricted 
to the functions (Ai, . . . A„) equals A^, as a consequence of (15. ip to (15. 4p : of course, 
this holds for the spectrum of every principal minor of the matrix B. D 



6. SDE FOR THE Dyson process on the spectra of two consecutive 

MINORS 

In this section we prove the formulas fl2.30p for the A- and //-SDE's, together 
with the generator (I2.32p . 



Proof of SDE li2.3(J\) and generator (I2.32p in Theorem \2.iA Using the same idea as 
in the proof of (I2.18P and (I2.20p in the Section |5l we choose, at time t, to conjugate 
the matrix B so as to have the form -Bbord of (I2.27P and let the matrix -Bbord evolve 
according to the Dyson process. We will consider only the first order effects on 
the A's and ignore second order effects. 

At first, we need to compute the (first order) variation of the A^'s as a function 
of the (first order) variation of the entries: 

Case 1: Consider the perturbation of -Bbord, using the notation (15. Sp for [dBij\, 
namely 



3 .18 I (6.1) Ebord + [dBij], for 1 < z < j < n - 1. 



Up to first order, one must compute the effect of the perturbation on each of the 
eigenvalues A^, by explicitly computing the characteristic polynomial of the bor- 
dered matrix (I2.27P with the extra non-diagonal perturbation; then, by neglecting 



3.19 
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the second order terms in dB, one findqj: 

(6.2) 



O=det(5bord+M5,,]-A/) 



Xl-^Xa+d\a 



n-1 



(6.3) =l[{f^e-X) 



'n-l 



(6.4) 



E 



^ X- ^i 



+ r-„ - A + 



2r,r 



(A-/ii)(A-/ij) 



b^E^a 



i;'(«.-p"' 



^=0 



AH-^Ac+dAa 



Setting A I— !■ Aq + (iAo, in this expression, shows that the product n"=i (/^^ ~ -^a ~ 
(iAo) is of the form (non-zero-function) + (function) x dX^^ whereas the second part 
gives, by Taylor expanding in A^, keeping in the expression first order terms only, 
evaluated by (13. 2p . and noticing that the 0*^-order term vanishes (again using 
dSSD), we find 



P'riK 



Pn-\\Xa) 



dXa + 



2rirj 



13-1 



(Aq, - fli){Xa - fJ'j) j^ 



J2dBl^{u.u;r = o. 



J^f'^ dB^f , with l<i<j <n-l, yields 



,(0) 
ij 



Finally, adding up the first order contributions from all the perturbations dB^^ + 

ij 



(6.5) 



dXr 



Pn-l{Xa] 

KM 



E 



2rirj 



13-1 



Y.^4 



'"'«.«•)'«. 



Case 2: For the perturbation [dBii\, with i = 1 
second order terms. 



. , n — 1, again neglecting the 



= det(5bord + [dBii] - XI) 

n-l 



Xh^Xa-\-dXa 

II (yU^ + 6eidBii ~ Xa- dXa) 



=1 



E 



■^ Xa + dXa — Hi — 6iidBi. 



+ (rn + SindBii) - Xa- dXc 



^Notice that 2RedB*.M,u* = dB*^u,u* + {dB*^u,u*y = 2J2iZo dB\f {u^u*)'^^\ using {ab)* = 
b*a* . Remember, for /3 = 4, quaternion multiplication does not commute and for the determinant 
of a matrix, use formula (I2.5p . 
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Upon expanding this expression as a function of Aq, jUj up to first order, noticing 
as before that the first part does not matter, and using again fl3.2p . this leads to 



(6.6) 
(6.7) 



PniK 



-dXa + 



-Pn-l(Ac 



dBii = for i = 1, . . . , n — 1, 

dXa + dB„n = for z = n. 



and thus summing up all the contributions coming from the dB^ for i = 1, . . . , n, 
one finds 



3.20| (6.8 



dXr 



Pn-l{\a) /^ r? 



E 



Case 3: For the perturbation [dBin], i = I, . . . , n — 1 



dBii + dB. 



M T^ ^-^nn 



= det(5bord + [dBin] - A/) 



Ah->Aq,+c/Aq: 



n-1 



n('"^-^) 



'n-\ 



^ rl + ri{uidBl^ + u*dBin)5ik _ 



,k=l 



>^- fJ'k 



Xh^\a+d\a 



Then, using UidB*^ + dBinU* = 2 J2e=o ""« (^^iL using formula (13.21) and finally 
summing up over all perturbations of the last row and column {1 < i < n — 1) 
yields 



n-l 



"3721] (6.9) 



dXr 



Pn-r{K) "y^ 2r,Y.tl4^dB\:! 



Then summing up the three contributions (16. 5p . (16. 8 p and (16. 9 p gives us the total 
first order contribution to dXa'- 



dXr 



Pn-iyXc 



2rirj 



y — 



/3-1 



^iuii^^ 



ydB;/[UiU, 



n— 1 2 n—\ 






> . 






Aq ~ ^i 
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We now set the SDE's (12. 8p for the dBa, dB-- into the equation obtained above, 
thus yielding, by (15. 5p . 



3:22] (6.10) dK = F^^\\)dt 



2 Pn-l\^a, 

1 Pk(K) 



E 






-^ {^a ~ f^i){^a — fJ'j) 



l<«<j<n-l r^/\ uc rj/ ^^Q 

n— 1 



4=1 ^ 



dbii + dh 



II I ^^nn 



+ ^Et4-E"."'* 



^ A„ -/i 



in 



£=0 



>, 



for some function Fq (A) to be determined later. Notice that in M, C and H, the 
norm |f | satisfies \vw\ = \v\.\w\ and |f | = |i;*|. Therefore, when \ui\ = 1, we also 
have \uiU*\ = 1, implying that 



0-1 



dbin := y^Ui dblJ and dbij := y ^{uiU*y^'db, 






1=0 



i=0 



are both standard Brownian motions on the sphere 5^^^; since they are differ- 
ent linear combinations, they are independent standard Brownian motions, and 
independent of dba, 1 < i < n. This is precisely formula (I2.30p of Theorem | 
namely 



sde 



(6.11) dX^ = Fi-\X)dt + .[l^^^ 



n-1 



y/2 TiTj dbij sr^ ^'^dbi. 



"" 'a<.l^n-i(^"-^*)(^"-^^) ' tr(^°-/^*)' ' tr ^--^' 



Ea/2 Tidbii, 
+ dbr, 



The SDE for the Dyson process induced on the {n — l) x (n — 1) upper-left minor 
is given by the first formula of Theorem 12.11 with n ^^ n — 1 and A H- /x, yielding 
the formula in (I2.30p . Therefore the product of the SDKs in (I2.30p . together with 
identity fl2:28|) yields 



3.23 



(6.12) 



uAjUyUj 2 J „_]^(Aj 



dt 



13 p^(A,) \\i-^lj 



2 1 Pn-l{\i)Pn{^J^j) 

'n\-Nfn{\)Pn-My 
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Moreover, 

= A^yXa^'^/) — Aq^^(A^) — A^^^(Aq,) 
= 2 (coefficient of — — — in A J 

OXaXi3 

-Is 

and similarly, 

[^ (6-13) ^r = /- 

These identities can also be computed from the expressions (16. lip of dXa in terms 
of the \i,fij, as done in the remark below. From Ito's formula l5. 5 [ it then follows 
that 



(ctAi, . . . , dXn, «/ii, . . . , dfln-l) 

= (^DysAl, . . . , ^DysAn, ^Dys/^l, • • • , ABysfJ'n-l)dt + 0"(A, ^)dbt, 

where, according to (I2.35P and (I2.2ip . 

^Dys(Ai) = AxiX^) = -Xi + y] :-, 

3.23"| (6.14) '^" ^ 

AuysifJ'i) = A^ifii) = -fti + y~] , 



estabhshing the form of #(A) in flOTj) . thus yielding flOOD . Identities flOOj) . 
(I6.12P and (l6J.3p . together with Ito's formula (15. 5p . then establish the formula 
fl233|) for Ax^ □ 



As an alternative way, (|2.35p and (I6.12p suffice to establish (|2.33p . with (|6.12p needed to 
establish the coupling A\^. 
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3.23' 



SDE-lambdaO ' 



Remark 6.1. Note that the identities fl6.13p can be computed as well from the SDE 
(16.111) . using residue calculations: 

(6.15) 

2 I Pn-l{^a)\ I Pn-l{^-y) 



di 



/3 V P^iK) 



^^(Ao 



2 rfr'^ 



vl<i<i<n-l (^Ci- ^-di^a- ^'j){^'y- l^i){\- (J'j) 
n— 1 A n— 1 ^ 2 



+E 



4 



^Ett^ 



2n^ 



^ (Aa-/ii)2(A^-/ii)2 -^ (AQ-/ii)(A^-/ii) 



/3 "^ 



dnidfjij 
dt 



—6ii. 
/3 '' 



We now turn to the proof of Corollaries 12.31 and 12.41 

Proof of Corollary \2.'J\ . Note, using logarithmic derivatives, that 

(invariant measure fl2.34p )/3^ (invariant measure (12.341) ) 



-1 



is a quadratic polynomial in /3, which by Theorem 12.21 vanishes for /3 = 1,2,4 
and thus it vanishes identically in /3. That the process restricted to A or /i is the 
standard Dyson process follows from the form of the generator AJ . 

U 

Proof of Corollary \2.4\ In order to study the stochastic behavior of /ic — A^ and 
Xa+i — fJ'a when fia gets close to Aq, or Aq,+i, one rewrites the Brownian part of 
d{Xa — /ia) as follows: 
(6.16) 



—Brownian part of d{Xa — fie 



-Pn-l(Aa) 



/ 



Ki>^ 



+ 



V 

Pn-l{X. 



l<i<j<n- 



\f2 TiTj dbij 



n-1 

+E 



r^dbr 



Pni^a) (Aq — fXa) 



^iXa-fii)iXa-fij) j^iXa-fiiY 
- 1 1 dbryr.. 



i^a 



Ea/2 Tidbin , ,, 
— ; h dbnn 

i=l 



Xa^ fii 



At first notice that for fi^ — Xa, one has, using the expression (12.281) for r 



fc' 



Pn-i{Xa) = 0{fia - Xa), T^ = 0{^/fia - Aq) and Tj = C(l) for i ^ a, 

from which one deduces that the first line on the right hand side of (16.161) has 
order 0{^\/ \ia — Aq). Using again (I2.28p . the second line of (I6.16P multiplied with 
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-Pn-l(/^a)-Pn('^a) 



p:^{K) (A. 



1 



/ioj 



hreeminors 



-Pn-l(Aa) — -Pn-ll/^a) \ ( Pn{l^a) — Pn{^ 
Aq /Iq 
0{Ha - A„), 



^Q, Ac 



^:-i(/i.)^:(A.) 



Then 



dt-part of (i(/iQ, — Aq) 



I /i.Q: — Xq 



c„_i ^^a l^j ,<,<„ Aq, Aj 






l<j<n 



A^CK — -^O 



E 



/ij - \j 



+ 



>o, 



i<j<„_i (Aq - ^^j){^a - Aj) A„ - Ac 

which follows from the inequalities (for 1 < j < "n. — 1), 

^^j ~ Aj > 0, (Aq, — fJ'j){^a — Aj) > and A„ — Aq, > 0. 
This proves the first relation (I2.36p . while the second one is done in a similar 



way. 



D 



7. The eigenvalues of three consecutive minors 



In this section we shall prove Theorem 12 .St which affirms that for the Dyson pro- 
cess the joint spectra of any three consecutive minors is not Markovian, although 
the Markovianess of the spectra holds for any one or any two consecutive minors. 

Note that given an Ito diffusion Xf G M", with stochastic differential equation 
dXf = a{Xt)dt + a{Xt)dbt, as in (15. ip . and generator A, the process restricted to 
Yi = ifi{X),l < i < i is not Markovian (at least for generic initial conditions) 
if the generator fails to preserve the field of functions J^{Y) generated by the 
{¥,,..., Ye) ■.= {viiX),...,MX)), i.e. 



4^ (7.1) 



AT{Y) ^-F(r). 



and provided the diffusion does not hit the F-boundary of the domain. 



Proof of Theorem \2.5[ In order to show the non-Markovianess of 

r := (A,/i,i/) = (spec 5, spec E^'^-^), spec S("-2)) 

it suffices to find a function, such that the function, obtained by applying the 
Dyson-generator to it, is not a function of (A,/i, z/). We pick a function of the 
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product form xy = g(r)h{r), where 



n-2 



X := g{T) := N^ Bu and y := h{T) := det B 



i=l 



672] (7.2) 



are two independent functions. Then, according to formula (15. 3p 

dxdy 



ABysXy 



dt 



+ xAoysy + yAjiysX. 



Since x and ^Dys^^ are functions of z/ only and since y and A^ysy are functions 
of A only, xAuysy + l/-4,Dysa^ is a function of (A, z/) only. Therefore, to establish 
non-Markovianess of (A,/i, i/), it suffices to show that -^ is not only a function 
of (A,/i, z/). Since, by Ito's formula (15. 5p . 



n-2 



(ix (iy = 2_, dBii 2_, ^^ — c^-B 



95 



jj 



jj 



^ gadet5^^(,) 






n-2 



2 , ^9det5 2^^ , , ,„,, , 

T^^^Z^ g^.. = ^ 2^det(mmorii(S))(it, 



/3 



it suffices to show that the right hand side is not a function of (A, /x, v) only. Here 
minorjj denotes removing row i and column i of the matrix. 

For example in the case /3 = 2, n = 3, this amounts to showing that the deter- 
minant of the lower-right 2x2 principal minor of B is not a function of (A, /x, v) 
only; to do this, it is convenient to reparametrize the matrix as 



B 



I B 

Pse 

\p2i 



11 



Pse 
B22 



^Vi 



-IV2 



oim 



pie 



'irix 



P2e 
Pie 

B33 ) 



Using the following formulae 

Bii 
B22 
B33 



^1, 

Ail +yU2 - Vi, 

Ai + A2 + A3 - /ii - P2 

(/i2-Z/l)(z/l - Pi), 

the lower-right 2x2 principal minor of B reads 

det(minorii(fi)) = 522^33 - p\ 

= {pi + P2- «^l)(Al + A2 + A3 - /Xi - P2) 



pI 



pI 



CONSECUTIVE MINORS FOR DYSON'S BROWNIAN MOTIONS 31 

One observes that 

= det5- A1A2A3 

3 
= BUB22B33 - A1A2A3 - ^ p. 5ii + 2pip2P3 cos(?7i + r]2 + r/3) 

= Fi{X,ij,,iy) - pIpi - pUpi + P2 -1^1) 



+ 2pip2 V(/i2 - z^i)(z^i - /ii) cos(?7i +r]2 + m) 
and 

= TrB^-iXl + Xl + Xl) = J2Bl + 2J2p'-i>^i + >^l + >^ 



I) 



= F2{X,p,iy) + 2{pl + pl), 

where Fj(A,/i, i/) are functions of the spectral data (A,yU, z/). Upon solving these 
two equations in pi and p2, one notices that, in particular, pi is a function of 
cos (771+772 +'73) and the spectral data (A, p, z/), hence showing that det(minorii(i?)) 
is not a function of (A,/i, z/) only; thus the same is true for Aj^ysXy. This proves 
that AnysXy does not belong to the field of functions depending on (A, p, z/). 

More generally, by a perturbation argument about i?("~^) = diag(/ii, . . . , Pn-i), 
one shows similarly that 

n-2 

y^ det(minorii(5)) ^ J'(A, /i, z/), 

i=l 

for /3 = 2 and 4. 

Finally, the boundary of the process (A,/i, z/) is given by the subvariety where 
some of the /ij's hit the A/s or the z^^'s; that is when Pn{pi) = or Pn-i{^k) = 
for some l<7<n — lor for some 1 < k < n — 2; r'^{p, u) = for 1 < j < n — 2. 
From Corollary 12.41 on sees that the process never reaches that boundary. This 
ends the proof of Theorem 12. 5[ D 
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